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ABSTRACT 

Exoplanets are often found with short periods or high eccentricities, and multiple- 
planet systems are often in resonance. They require dynamical theories that describe more 
extreme motions than those of the relatively placid planetary orbits of the Solar System. 
We describe the most important dynamical processes in fully-formed planetary systems 
and how they are modeled. Such methods have been applied to detect the evolution of 
exoplanet orbits in action and to infer dramatic histories from the dynamical properties of 
planetary systems. 



1. INTRODUCTION 

After a planet has formed via giant impacts and has 
outlasted migration torques from the gaseous disk, per- 
ils still await. Dynamical instabilities among the plan- 
ets of the system, long-term orbital changes due to a 
companion star, and tidal interactions with the host star 
could all eject the planet or toss it into the host star 
None of these would threaten planets if they contin- 
ued on Keplerian orbits (as described in chapter 1 of 
this volume), in which the planet and star are consid- 
ered as point masses, orbiting each other according to 
Newton's approximate theory of gravity, in isolation 
from other bodies. To model these more interesting 
interactions, we explore non-Keplerian dynamics, an 
extension of orbital theory to (1) the astrophysical two- 
body problem in more detail, including non-spherical 
bodies and relativity theory, and (2) true planetary sys- 
tems: systems with more than just one planet and one 
star. 

Although most dynamical concepts were originally 
designed to describe the Solar System, the discovery 
of exoplanets has channeled research in new direc- 
tions. Previously most analytical and even numerical 
techniques required certain quantities (mass ratios, ec- 
centricities, mutual inclinations) to be small. These 
approximations are appropriate for applications in the 
Solar System, but more general methods are needed to 
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describe the bewildering variety of exoplanetary sys- 
tems. Here we describe these methods: both new di- 
rections of analytic theory and an introduction to nu- 
merical work. 

Orbital motion departs from fixed Keplerian el- 
Upses on a wide range of timescales, from orbital pe- 
riods to stellar lifetimes. These variations have a wide 
range of magnitudes, from slight deflections to com- 
plete reorientations and ejections. Slight deflections, 
due to planets passing one another, may be too small 
to detect, apart from sensitive transit timing measure- 
ments. Secular interactions and dissipative effects may 
be too slow to detect in individual systems, but they 
may be infeiTed by statistical studies of populations. 
Ejections may be too rare to see in action, but they 
are expected from numerical modeling. Resonant in- 
teractions, however, can produce substantially non- 
Keplerian motion over a timescale suitable for obser- 
vation, and resonant systems can stably persist for the 
star's entire lifetime. 

The plan of this chapter is as follows. In |j2] we be- 
gin with equations of motion for planets that include 
tidal distortion, rotational oblateness, effects of relativ- 
ity, and gravitational interactions among planets. The 
framework for solving these equations is presented, in- 
cluding introductions to both coordinate systems and 
numerical algorithms. In |j3j the phenomena which 
arise from these equations are discussed. We empha- 
size phenomena which are important for known sys- 
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terns, as well as types of orbits that could exist but 
are not yet observed. Next, in [JH we turn to obser- 
vational highlights which have used and continue to 
challenge these concepts, showing how theories and 
systems in nature have enjoyed a symbiotic relation- 
ship. Finally, ijS] discusses observations that would be 
particularly useful for constraining theories of dynami- 
cal evolution. With this structure, some individual top- 
ics in dynamics reappear in several different sections; 
so we provide Table[T]as a topical guide. 

2. EQUATIONS OF MOTION AND NUMERI- 
CAL METHODS 

2.1. Astrophysical Two-body Problem 

First, let us examine the astrophysical two-body 
problem including non-Keplerian effects. The masses 
of the star and planet are TO,t and nip, respectively, and 
the orbital elements refer to the displacement vector r 
(of magnitude r) of the planet relative to the star, as in 
Chapter 1 of this volume. The equation of motion is: 



r 
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(1) 



where / is a force other than that of mutual gravity 
of point masses, so f = yields Keplerian motion 
(Chapter 1, equation 4). 

2.1.1. Relativistic effects 

The lowest-order post-Newtonian effects may be 
implemented using the force: 
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MardUng & Linll2002l) for f in equation ([T). We will 
examine one dynamical consequence of this force, ap- 
sidal motion, in i jS.l.ll Alternatively, a potential that 
mimics lowest-order relativistic effects, which is es- 
pecially suitable for analysi s within a Hamiltonian 
framework, was given by (ISaha & Tremaine 1 19921 
eq.31). 

2.1.2. Effects of non-spherical bodies 

The orbital parameters change when we consider 
the bodies not as point masses, but as physical objects 



capable of distortion and internal energy dissipation. 
Tidal effects on the orbit become more and more pro- 
nounced as two gravitating bodies of finite extent get 
closer to one another. Tides have apparently caused 
many exoplanet orbits to become circular within about 
0.1 AU of their main-sequence stars. 

The tidal force of each body distorts the poten- 
tial energy surfaces of its companion. For stars and 
gaseous planets in hydrostatic equilibrium, the sur- 
faces of constant density will settle to these new 
equipotential surfaces, and the distortion of the body 
itself modifies them further, until a self-consistent so- 
lution is obtained. For solid planets, the rigidity of the 
material also comes into play. These properties deter- 
mine the planet's tidal defo rmability an d are summa- 
rized by a Love number fc^ ( Lovel[l9Tl ). which is the 
amplitude ratio of the quadrupolar potential due to the 
deformed body to the tidal potential imposed on the 
body, evaluated at the surface of the body; a table of 
typical values of may be found in IMardling & Lin 
(2004). Below we shall also denote the Love number 
of the star as fci_,t, which is twice the apsidal motion 
constant, the conventional parameter in the literature 
regarding eclipsing binaries (Ster ne 1939). 

First, the star raises a tidal bulge of the planet with 
a size oc r^"^. This bulge creates its own external field 
which falls off like r^^. Acting back on the star, these 
radial scalings combine to augment the radial gravita- 
tional force (per unit mass) with a term having steep 
radial dependence: 



fx = 



(3) 



Second, the rotation of the star causes it to become 
oblate. Its degree of oblateness depends on the square 
of the stellar angular rotation rate il* divided by its 
break-up angular rate, ^/GmT/B^. The quadrupolar 
potential outside of the star scales as r^^ . Thus the 
figure of the rotating star induces an extra force 



(4) 



The forces of the foregoing two subsections may 
be included in numerical integrations to simulate more 
realistic orbital behavior than Keplerian ellipses. They 
may also be used in analytic calculations to determine 
orbit changes on a longer timescale, results of which 
we quote in i]3.1l 
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2.2. N-body equations and coordinates 

When a star hosts N (> 1) planets, gravitational 
interactions among the planets can affect their orbits 
in complex ways. Neglecting the effects discussed in 
^2.11 the equation of motion of planet planet i (of mass 
nii) is: 



- r, 



-12 

' 

(5) 

where each of the coordinates is referred to the cen- 
tral star, of mass toq = m*. The interaction terms 
in the sum over each of the other planets are the di- 
rect gravitational force (first term) and the indirect ef- 
fective force due to the bodies causing the star, and 
thus the reference frame, to accelerate (second term). 
In numerical work, these N second-order differential 
equations are most often transformed into a system of 
2N first-order differential equations in the quantities 
Tj and r^. In response to the motions of the planets, the 
star's position and velocity with respect to the barycen- 
ter of the system (the center of mass) are: 



Ro = -(^^mirj^/(^^TOi^ (6) 

i=l 1=0 

m,ri'j /(^Y^ m,) , (7) 



which are useful for self-consistent fits of data (see 

The coordinates of equation (|5]l are called astrocen- 
tric, because they refer each planet to the (moving) po- 
sition of the star. There are several other possible co- 
ordinate systems — barycentric, Jacobian, or Poincare 
coordinates — which are defined as follows. Barycen- 
tric coordinates refer all positions and velocities to the 
center of mass, and the only force is Newton's attrac- 
tive force between bodies (no indirect term). Jacobian 
coordinates are hierarchical, in which the positions and 
velocities of each of the planets is referred to the cen- 
ter of mass of all interior bodies. Poincare coordinates 



I Laskar & RobuteT^lQQS; also called democratic helio- 
centric coor dinates (Duncan etal . 1998) or canonical 
heliocentric ( Sussman & WisdomlboOll pg. 435)], on 
the other hand, take the positions of the planets to be 
astrocentric, but take the velocities of the planets to be 
barycentric. 

The various coordinate systems have different 



strengths and weaknesses, in both numerical integra- 
tion and in analytical studies, as follows. 

Astrocentric coordinates have a few drawbacks. In 
numerical integrations, all planets must be integrated 
with timesteps smaller than the innermost planet's pe- 
riod, because of the high-frequency forcing it intro- 
duces. Also, for theoretical studies in a Hamilto- 
nian framework, the astrocentric coordinates are rather 
cumbersome (Beauge et al. 20 07). 

Although barycentric coordinates are very simple, 
they do not take advantage of the fact that the cen- 
ter of mass is invariant. Therefore the position of the 
star must either be numerically integrated or updated 
at each timestep according to equation A reduc- 
tion by those 3 degrees of freedom is achieved using 
Jacobian or Poincare coordinates. 

iLee & Pealei (120031) have advocated using Jacobian 
coordinates, instead of astrocentric coordinates, when 
turning the solutions of multiple Keplerian radial ve- 
locity fits into a self-consistent N-body realization 
( m.Ai . Particularly for systems that are hierarchi- 
cal, these are the coordinates in which the planets 
perturb each other minimally on orbital timescales, 
so the independent-Keplerian model is satisfied best. 
Once an integration is set up in Jacobian coordinates, 
one might want to do the integration in astrocentric 
coordinates, so that, e.g., transit times of the outer 
planet can be easily calculated ( ^12. 51 ). The hierar- 
chical structure of Jacobian coordinates is not always 
physically well-motivated, e.g., for systems with Tro- 
jan orbits. Most authors take the relevant mass bind- 
ing the j*^ planet of mass mj to the mass interior to 
it TOintcriorj to be TTij + mi„tcrior,j (in generalization 
of the mass term eq.[T]), but in symplectic integrations 
(described in m.3i it can be more efficient to use M, = 

TO^^interior. i /"^interior. i - 1 dWisdom & Holmanll991 

Murrav & Dermottlll999h . 



Poincare coordinates have positions and velocities 
with different origins, so they are not conceptually 
based on elliptical motion. However, analytic studies 
can be cleaner, as the resulting Hamiltonian has rather 
clear symmetry. 

For each coordinate system, each planet has "oscu- 
lating" orbital elements which correspond to the orbit 
the planet would have if all the interaction forces are 
ignored. For instance, orbital elements based on the 
astrocentric coordinates have a simple physical inter- 
pretation: if {N — 1) planets of the system suddenly 
disappeared, a single planet would be left orbiting the 
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star with those orbital elements. Osculating elements 
have a draw-back, which is that the physical, observed 
properties, which include the continual perturbations 
of the other bodies, do not correspond to the osculat- 
ing value. For instance, the orbital period observed 
with transits or radial velocities is not the osculating 
period in an astrocentric coordinate s ystem (for a con- 
Jperraz-Mello et al.ll2005h . 



version, see 



For explicit forms of t hese coordina t e syst ems, 
and a fuller discu ssion, see iBeauge et al.1 (120071) and 



Morbidellil (120021) 



2.3. Numerical integration techniques 

Here we summarize some of the numerical tech- 
niques and codes that are used in exoplanet dynamics. 

Currently, the fastest algorithms that are reliable 
for long-term studies are symplectic integrators, in 
which separate integrations delimiting a certain vol- 
ume of phase space continue to delimit the same vol- 
ume as the integration proceeds. This property auto- 
matically respects the Hamiltonian nature of the grav- 
itational pro blem; e.g., energy is not sec ularly lost or 
gained. The IWisdom & HolmanI (1199 lb integrator is 
the most well-known of the symplectic integrators. It 
has been supplemented with algorithms that can han- 
dle close approaches and collisions, which forms the 
backbone of the publicly- available codes Mercur^F] 
jLevison & DuncanI [l994 ) 



backbone or the publicly-i 
Jchamberl lT 999l). Swim I 



and SvMBA^ dOuncan et al.lll998h . A symplectic al- 
gorithm was also built specifically to h andle stellar 
mass companions ( [Chambers et al. [l2002h. either with 
planets orbiting one star of the pair (satellite-type) or 
both stars (planetary-type). The speed of symplectic 
integrators allowed some of the first long-term inte- 
grations of all the planets of the solar system, show- 
ing th eir orbits to be chaotic (e.g..l Sussman & Wisdom 
19921) . High-order symplectic schemes have been 
used to follow the whole S olar System for its lifetime 
dLaskar & Gastineaull2009h . 

There are several frequent ly-used, but non -symplectic, 
algorithms. RADAU (15) bv lEverhartl(ll985h is a very 
high-order method of integration by Radau quadra- 
ture — a prescription for times to evaluate forces 
and weights to apply to them — and it is included 
in Mercury. Ge neral-purpose inte grators from Nu- 
merical Recipes dPress et al.lll992l e.g., the popular 
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Burlisch-Stoer) and the GNU Scientific Library [e.g.. 
Embedded Runge-Kutt a Prince -Dormand (8/9 order), 
as advocated by Frege auet al.l 120041 have also been 
set to work on problems in exoplanet dynamics. These 
integrators are particularly useful for short-term, high 
precision work. 

If mean motion resonances are unimportant, one 
can derive approximate secular equations of mo- 
tion by analytically averaging each planets' grav- 
itational effect over its orbit. This procedure re- 
sults in differential equations for the orbital ele- 
ments, which can be numerically integrated much 
faster than the positions themselves. iMardling & Lin 
(I2OO2I) have presented equations to evolve three- 
planet systems, as well as separate equations for 
the system once the inner planet's or bit is averaged. 
'Eggleton & Kisele va-EggletonI (I2OOII) have presented 
equations for the secular problem for three-body sys- 
tems where the semi-major axis ratio is large, and 
the outermost orbit dominates the angular momen- 
tum budget; this approach is applicable to planets in 
binaries. Higher order expansions in the ratio of semi- 
major axes, with an emphasis on planetary systems, 
have a lso been derived (Ford et al. 2 000; Lee & PeaI3 
20031 iMigaszewski & Gozdziewskil boosi) . Taking 



no restrictions on relative size of orbits or orbital 
elements, one can compute by brute force the time- 
averaged torque acting among a set of nearl y-Keplerian 
orbits, which is called Gauss's method (e.g.. lTouma et al 



20091) . 



Finally, a general consideration, which is surprising 
at first, is that codes are usually much slower at print- 
ing data to files than computing data. For instance, 
if the position of the planets is printed once per orbit, 
the print-outs will generally completely dominate the 
runtime. It is generally true that performance can be 
limited by the weakest link (the slowest part of the al- 
gorithm), so one must think about the program as a 
whole, including the reads and writes, when trying to 
speed up an algorithm. 

2.4. N-body fits to data 

To directly detect non-Keplerian motion, a frame- 
work must be constructed for comparing the outputs 
of numerical algorithms to data. We describe the con- 
ventional approach here. 

The astrocentric coordinates for planets are = 
{xi,yi, Zi); the barycentric coordinates of the planets 
are Ri = (X^, 1^, Zi), and the barycentric coordinates 
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of the host star are Rq — (Xq^Yq^ Zq), which obey 
equation (|7]i. The standard coordinate system (e.g., 
Chapter 1) has the Z-axis pointed away from the ob- 
server. Let X — Y plane be the plane of the sky, which 
passes through the barycenter Then aligning the X- 
axis with the North direction and the Y-axis with the 
East direction forms a right-handed coordinate system 
with Z. The ascending node ft is measured East from 
North, and the argument of periastron u is measured 
from the sky plane. These coordinates are useful for 
fitting data from essentially all techniques (neglecting 
the rectilinear motion of the system as a whole); 

1. radial velocity of the host = —Zq, 

2. astrometric position of the host = {Xq, Yq), 

3. transits and secondary eclipses may be modeled 
using the light curve equations of Mandel & Agoll 
(l2002h . setting the projected distance between 

^ (see 



the centers, called d, equal to ^ya 
also i l2.5l below). 



4. direct images of planets have a sky-offset from 
their host of {xi,yi), 

5. the times of arrival (TOA) of pulses of pulsars, 
or phases of pulsating stars, or eclipses of close 
binaries, that are orbited by planets are delayed 
by 6t = —Zq/c, etc. 



Therefore a common coordinate system can be used 
to fit different types of datasets simultaneously, as has 
already been done for radial velocity plus astrometry 
feean & Seifahrt 2009; Wright & Howard 2009) an d 
for radial velocity plus transits ("e.g. JWinn et al.l 2005 ). 

The values of G and tuq are not known to nearly 
the precision of their combination. A result is that 
dynamical fits reported in the literature are typi- 
cally not reproducible if the masses are reported in 
physical units (e.g., grams). One solution is to re- 
port all masses in units of solar masses. Alterna- 
tively, one could abide by the International Astronom- 
ical Union's conventional value of fc = Gtuq = 
0.01720209895 AU3/2 day-^. 

There are numerous methods for converging on the 
solution and evaluating uncertainties. Markov Chain 
Monte Carlo (Ford 2005, 2006), Levenberg-Marquard t 



dLaughlin & Chambersl20oTllRivera & Lissauej200"ll) 



tempering; lGregorvll2007h have all been used exten- 
sively, and these references describe the implementa- 
tions and advantages of each method. 

The magnitude and/or timescale of perturbations 
scale as powers of the mass ratios, so the true planetary 
masses in systems with measurable non-Keplerian mo- 
tion are accessible (i i4. 114.21 1. Moreover, often only a 
subset of orbital parameters which fit the observations 
actually result in long-term stability, so additional in- 
formation about the system is accessible by requir- 
ing long-term stability. This idea was pursued for the 
first multiplane t system by Riv era & Lissau er (2000|); 
Stepinski et al.l (l2000l) . and it ha s been particularly 



useful for resonant systems; 
(l2008bl) for an overview. 



see 



Gozdziewski et al 



genetic algorithm d Gozdziewski et all 2005 ), and mul 
tiplexed simulated anneahng (also known as parallel 



2.5. A focus on transit variations 

Now let us focus on dynamical fits to transit data, 
which have the potential to discover small planets and 
to reveal the detailed dynamical properties of multi- 
planet systems (see ij5]l. The data begin as flux mea- 
surements of a star as a function of time; the star is 
dimmed as the transiting planet covers some of the 
stellar surface (see chapter 4). Encoded in this time 
series is information about the position of the planet 
relative to the star, so integrators based on astrocen- 
tric coordinates are the most natural for this applica- 
tion. The data allow for precise positions and times to 
be observed, so high precision in coordinates and ar- 
bitrary timesteps is required. Typically, the flux mea- 
surements are used to derive transit parameters — mid- 
transit times and durations (or impact parameters or 
inclinations) — as a function of transit number, which 
may be compared directly to a numerical integration. 

Theoretical mid-transit and contact times for the i*'' 
planet may be computed by integrating forward and 
backward in time by the Newton-Raphson method, 
seeking roots of functions of r^.; = {xi, yi), the rel- 
ative separation vector of the planet and star on the 
plane of the sky; see Figure [T| This method is de- 
scribed next. The only requirement for the following 
algorithm is that the sky-projected trajectory has a ra- 
dius of curvature larger than i?^ + Rp, which is easily 
fulfilled in practice. 

For each transit the mid-time is found first, as fol- 
lows. The first step is to advance the integrator in time 
to the vicinity of an observed transit, for instance, by 
stopping an integration once an observed transit time 
has been passed, or advancing the appropriate num- 
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ber of nominal orbital periods from a previous transit. 
Mid-transit times may be found by minimizing Ir^.^j. 
Taking the derivative of this quantity (expressed in Xi 
and yi) and setting it to zero, we find that minimizing 
|rs i| amounts to solving 



XiXi + yiiii = 0. (8) 



The first guess at the time of transit must be within 
about 1 /8 of a period, or the algorithm may converge 
to an adjacent local maximum of |rs,i|. The solution 
of equation[8]is found by 5-10 iterations of moving the 
integrator by 



8t 



fdg 



where 



— = + x^Xi +yi + y^y^ 
ot 



(9) 



(10) 



according to equation (|8]l. Once 6t is below the re- 
quired accuracy, the mid-time tniid and position r""*^ 
of that transit have been found. This method is compu- 
tationally considerably faster than searching for a min- 
imum of |rs i| directly, as root-finding is a simpler op- 
eration than minimum-finding. 

Next, before solving for times of contact, we can 
determine if they exist for this particular transit. (See 
chapter 4 for the definitions of points of contact and 
for the numbering scheme.) If |r""^| < R^, — Rp, 
then all four times of contact exist; the trajectory is not 
grazing. If R^ - Rp < \rf^f^\ < R^, + Rp, then the 
trajectory is grazing and only first and fourth contact 
exist. If |r™^| > _R* + Rp, the planet is not transit- 
ing at all and searching for times of contact is a waste. 
To search for times of contact, first advance the inte- 
grator to tniid T where — is taken for first 
and second contact, and + is taken for third and fourth 
contact. Next, solve the equation 



h{xi,yi 



- (R, ± Rpf = 0, 



(11) 



where + is taken for first and fourth contact, and — is 
taken for second and third contact. As before, this is 
done by iteratively driving the integrator by 



where now 



St 

dh 
'di 



-h 



\dt) ' 



(12) 



(13) 



It is important not to start the search near mid-transit, 

dh 
dt 



Otherwise ^ = 2g and equation (fTZt will ask 



for an enormous jump: preemptively advancing to near 
or beyond the times of contact (as above) avoids that 
fate. As usual with Newton-Raphson's method, sensi- 
ble bounds should be placed on the search for the root. 

After times of contact and mid-transit are com- 
puted, they need to be corrected for the finite travel 
time of light. Thus the apparent time of an event de- 
pends on the distance between the bodies and the cen- 
ter o f mass of the system, along the line of sig ht (see. 



e.g., 
pose 



Knutson et al. 2007: Bakos et al. 2009; 



Sup- 



for instance, that the ephemeris of mid-transit 
times for a transiting planet is well-established, and 
the values of eccentricity e and longitude of periastron 
w are known. Then the time of secondary eclipse could 
be predicted, but it will actually be observed a time At 
later due to the finite speed of light, where: 



At^ 



2a mi 



1 



c (m^ + mp)2 1 - sin uj 



(14) 



This equation is derived by delaying the image of each 
body according to its line-of-sight distance, then solv- 
ing for the time at which these images cross, i.e., the 
mid-time of secondary eclipse as observed. This delay 
can be large for transiting planets, e.g. At w 160 s 
for HD80606b, a systematic effect that is comparable 
in magnitude to the measurement error of ^ 260 s fo r 
an individual secondary eclipse (Laug hlin et alJE009l) . 
The reason this effect has an observable magnitude is 
because it originates within special relativity, and so 
it scales as 1 /c. Effects of general relativity, such as 
those embodied in equation (O or the effe ct of cur- 
vature on light propagation ( Shapiro! 19641) . scale as 
1 /c^, and are thus harder to detect. 

A final step is needed to compare these theoretical 
event times to the observed data. The most important 
effect of a second body on transit times is to change 
the true period (as mentioned in m.2i . so in practice 
the times are scaled to match the observed period and 
shifted to match the observed epoch of transit. Then 
variations on top of that simple linear ephemeris can 
be used, e.g., to detect or set limits on the presence 
of perturbing bodies ( il5.ll ). Care needs to be taken if 
radial velocities are being simultaneously fit, because 
such data specify a particular velocity of the star cor- 
responding to a particular time, whereas transit data 
specify a particular time corresponding to a particular 
relative position of the planet to the star. 
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Fig. 1. — Diagram of a transiting planet at the contact 
and mid-transit points of its orbit. 

3. DYNAMICAL PHENOMENA 

3.1. Astrophysical Two-Body Problem 

3.1.1. Periastron advance 

In m.l\ we encountered several extra forces which 
modify realistic two-body motion from Keplerian el- 
lipses. Let us take a perturbation theory approach, in 
which the force is calculated over a Keplerian o rbit, to 
find how that orbit itself evolves ( lBurnslll976l) . Now 
that the effective force of gravity is no longer of the 
form l/r^, the ellipse does not close, and the perias- 
tron advances, which amounts to a reorientation of the 
orbit within its own plane. We shall calculate this pre- 
cession rate due to the perturbing forces of relativity, 
tides, and rotational oblateness. 

Periastron advance is the relativistic effect which 
changes orbits on the shortest timescale, and averaging 
equation ^ over an orbit, we determine its angular 
rate to be: 



WGR 



05/2^2(1 _g2) ■ 



(15) 



This effect causes an additional 43 seconds of arc per 
century of precession for Mercury (in addition to the 
precession caused by the other planets), which was the 
famous first hint that nature obeyed Einstein's equa- 
tions. Higher order corrections to this precession rate, 
precession due to the star's spin (the Lenz-Thirring ef- 
fect), orbital decay due to gravitational wave emission. 



and other relativistic effects are all negligible for exo- 
planets. 

The force due to tidal distortion of the planet 
(eq. ||3|) causes apsidal motion at the rate: 
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where n = 27r/P is the mean motion ( Sterne 19391) . 
This effect is generally much bigger than that of the 
tide raised on the star by the planet, and for hot Jupiters 
with periods less than thre e days, it can dominate al l 
other precessional effects (Ragozzine & Woljl2009l) . 
For physically smaller planets, and for Jupiter-type 
planets with periods > 3 days, relativistic precession 
(eq. ifTSl ) typically dominates. 

Rotational distortion gives rise to a force (eq. HJ) 
that causes a periastron advance rate: 



n/CL,* 1 + TOp/TO* ^ 



(1 — e^)^ \ n / \ a 



(17) 



Around fast-rotating and large stars (i.e., young or 
early-type), this effect can dominate the others. If the 
stellar spin is misaligned with the orbit by an angle 
equation ( fTTI i requires an extra term (5 cos^ ip — l)/4; 
for spin-orbit angles satisfying 63.4° < ip < 116.6°, 
the apsidal motion is retrograde. With spin-orbit mis- 
alignment, the nodal angle also precesses; equations 
for the coupled spin and orbital motion are given by 
Eggleton & Kiseleva-Eggleton (2001). 

Of course, the star has a tidal bulge and the planet 
has a rotational bulge as well, but these never con- 
tribute substantially to the total precession. 

3.1.2. Tidal dissipation 

Tidal energy is converted to heat when a tidal bulge 
rotates through a body or varies in amplitude, due to 
the material's resistance to shearing motion. First, the 
dissipative torque changes the rotation of the planet to 
a rate at which the time average of that torque vanishes. 
At this spin rate the time average of the shear, and the 
energy dissipation rate, is minimized. In a fixed, circu- 
lar orbit, the spin angular velocity equals the orbital an- 
gular velocity and the obliquity is zero, so in the frame 
corotating with the perturber, the tide is no longer time 
variable, stopping energy loss. In an eccentric orbit, 
the spin will either se t tle at a pseudo-synchronous stat e 



(lPeale& Gold! 1 19651: lHutill98U iLevrard et al.1 120071) . 



or be trapped in a spin-orbit resonance (of which Mer- 
cury is the prototype); the latter is only possible if the 
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body has a permanent quadrupole moment due to its 
rigidity, and is therefore not expected for gas giants. 
For rocky planets with dynamically-important atmo- 
spheres (of which Venus is the prototype), the picture 
can be qualitativel y different, inc luding up to four sta- 
ble rotation states ( Correia et al ]|2008). 

On a longer timescale, the eccentricity damps. The 
correlation between eccentricity and orbital distance 
(or period) is the main constraint on tidal theory for ex- 
oplanets (see chapter 9). This damping can in principle 
be due to either dissipation in the planet or the star If 
dissipation in the star is important, eventua lly the plan 



etary orbit will decay into the star (e.g.. iRasioetaL 
19961: IJackson et al.ll2009t barker & Ogilviell2009l ). If 



dissipation in the planet is important, then it may have 
ingested more tidal energy than its own binding en- 
ergy. I n that case, gas giants could inflate or even 
disrupt (iGu et a l. 2003), and such heating on terres- 
trial plan ets would have significant geophysical conse- 
quences d WisdomlbOOSh . 

The physical causes of tidal damping for giant plan- 
ets are still poorly known; as yet no first-principle the- 
ory is efficient enough to damp the eccentricity of hot 
Jupiters, or to generate the inferred histories of satellite 



system s aro und the four So l ar Sys tem giants. iLin et al. 



(I2OOOI) and lOgilvie & LinI (|2004 discuss these mat- 
ters and review the literature. In the absence of such a 
theory, a phenomenological app roach has gained cur- 
rency jOoldreich & Soteij[l96 6^. A fraction l/Q of 
the tidal energy is dissipated per tidal forcing cycle (or 
per orbit, depending on the author). This allows dif- 
ferential equations for tidal damping to be derived, in 
which damping times scale w ith Q (Mardling & Liril 
2OO2I: iMatsumura et al.ll2008h . Empirical constraints 



on Q for close-in ga s giants have been worked out (IWu 
2OO3I: IJackson et al.i,2008; .Matsumura et al.ll2008h . 



3.1.3. Miscellaneous orbital evolution 

There are numerous other effects that can modify a 
planet's orbit about its star Here we simply list some 
of these effects, referring the reader to work that de- 
scribes them in detail. 

Close in to the star, the planet may be tidally 
stripped of mass. As the mass leaves the planet, it 
applies a torque on its orbit. The reaction of the orbit 
has been calculated for circular or bits as the plan et fin- 
ishes migration due to a gas disk dTriUing et al .|[T998i 
for moderate eccentricities a s the planet tida lly circu- 
larizes and perhaps inflates ( Gu et al. 20031) . and for 



eccentricities near 1 when the planet is shot near the 
star by eithe r a dynamical instability or a chance flyby 
dFaber etal.l Eo05). 

Once close to the star, the planet's atmosphere ab- 
sorbs and reradiates photons in preferential directions, 
which can lead to at most a 5% change in semi-major 
axis — enough to influe nce resonant co nfigurations 
with more distant planets ('Fabrv ckyl2608l) . 

A planet or planets may scatter and eject a sea of 
small bodies (planetesimals) after the main forma- 
tion phase, which leads to planetary migration. This 
effect was first wo rked out for the giant planets of 
the Solar Sy stem dFernandez & Ipl 1 19841: iMalhotra 



1993aL 119951). and has sinc e been apphed to exo- 



planets dMurrav et al.l 119981: iMorbideUi et al.1 12007 



Thommes et al.ll2008l) . For more on migration, par- 
ticularly in a gas disk, see chapter 14. 

Finally, far from the host star, pas sing stars may 
pertur b planetary orbits (see, 
2009h. 



e.g., ISpurzem et al 



3.2. Short-period, secular, and resonant interac- 
tions 

The interaction terms in the equations of mo- 
tion (eq. Q) lead to all the interesting behavior of 
N-planet systems, and here we show heuristically how 
short-period, secular, and resonant behaviors arise. 
(For a traditional expansion in terms of orbital ele- 
ments, using the so-call ed disturbing function, see 
Murrav & Dermottl[l999[ ^ Figure |2] shows how the 
Jacobian orbital elements of two planets, initially on 
circular orbits, evolve as the planets move from op- 
position (being on opposite sides of the star), through 
conjunction (lined up with the star on the same side), 
then back to opposition. As the planets approach each 
other, the mutual gravity moves them onto slightly dif- 
ferent orbits. At conjunction, the inner planet has been 
torqued forward, to a more distant, slower orbit; the 
outer planet has been torqued backward, to a closer-in, 
faster orbit. As the planets move through and recede 
from conjunction, these orbit changes are mostly re- 
versed. 

Now let us introduce moderate eccentricity to the 
orbits, and follow the system for many conjunctions. 
Due to the eccentricity, the paths of the planets at var- 
ious conjunctions are either converging or diverging, 
and the changes in orbital elements on either side of 
conjunction do not cancel as completely. At a single 
conjunction, this causes the orbits to transfer energy 
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Fig. 2. — Orbital element changes on the timescale 
of conjunctions, in a hypothetical system. The stellar 
mass is m©, both planets have mass IG^'^mQ, their or- 
bits are coplanar and initially circular, and they start on 
opposite sides of the star. Here and elsewhere, plan- 
ets are numbered by increasing semi-major axis, top 
panel: the mean longitudes of each planet (the two 
lines cross at conjunction), second and third panels: 
the planets' semi-major axes, which vary symmetri- 
cally about the conjunction,/oMrf/! and fifth panels: the 
planetary eccentricites, which receive a small kick. 



(the semi-major axes change) and angular momentum 
(the eccentricities and orbit orientations change). After 
multiple conjunctions, the behavior of the system de- 
pends on whether the periods are near a ratio of small 
integers. 

First, consider Figure |3] which shows a hypothet- 
ical system with a period ratio far from any ratio of 
small integers. Because of this property, the conjunc- 
tions sample all parts of both orbits rather equally. 
The semi-major axes exhibit no long-term changes, 
which means the energy of each orbit is conserved. 
However, the angular momentum of each orbit is ex- 
changed on long timescales, resulting in eccentricity 
variations. One way of seeing why this happens is by 
considering the time average of a planet over its or- 
bit, so that its gravitational effect is that of an elliptical 
wire weighted inversely to the Keplerian velocity at 
each position. Each of the planets respond to the other 
planets as if they were such rings. Because the poten- 
tial from such a ring is not time dependent, it produces 
a conservative force, and no energy can be exchanged: 
semi-major axes and periods may not change. How- 
ever, the lopsided rings do torque each other, and this 
corresponds to angular momentum (and thus eccen- 
tricity and inclination) changes. For instance, two ec- 
centric, coplanar planets will undergo periodic oscilla- 
tions in eccentricity which are 180° out of phase from 
each other. For more on this topic, called called secular 
evolution, see [J33]and ^4.3\ 

Next, consider Figure |4] which is a hypothetical 
system with periods very close to a ratio of small in- 
tegers (2.01:1). In this situation, called a mean-motion 
resonance, conjunctions occur at the same part of the 
orbit many times in a row, and the change in orbital 
elements builds. One may consider a changing period 
the hallmark of a mean-motion resonance. Along with 
period changes come eccentricity oscillations, which 
can be rather large over only tens of orbits. For non- 
coplanar planets, the distance between the location of 
conjunctions and the intersection of the orbital planes 
affects their dynamics, so resonances can also involve 
inclinations and not only eccentricities. In general, the 
angles that dictate the behavior of the resonance are 
called critical angles, and they have the form: 

(f> = jlAi+j2A2+j3^^1+i4^572+j5^^1+i6^^2, (18) 

where each planet has a mean motion of A, a longitude 
of ascending node ft, and a longitude of periastron to. 
The j values are integers obeying J^ji = (called 
the d'Alembert relation), which is required by the in- 
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variance of the system's behavior to the arbitrary ref- 
erence direction from which angles are measured. At 
the very center of each resonance, where the planets 
come to conjunction at exactly the same point in their 
orbits, the periods are constant and precession rate is 
constant. In this case, there exists a slowly rotating 
frame in which the motion of each planet is perfectly 
periodic, yet not perfectly elliptical. 

Individual resonances can help keep a system sta- 
ble. For instance, when the critical argument for the 
interior 2:1 resonance has zero libration amplitude 
(9 — 2A2 — Ai— n7i = 0, ^ = 0), we may rearrange the 
equation to read Ai — A2 = A2 — rui, which shows that 
when the two planets are at conjunction (Ai — A2 = 0), 
then the inner body is also at pericenter (Xi~wi — 0), 
so a close-approach is avoided. Conversely, whenever 
the outer body is at the azimuthal location of the inner 
body's apocenter (A2 — tui = tt), the two bodies are 
farthest apart (A2 — Ai = tt). This argument, applied 
generally to mean-motions resonances, is called a res- 
onance protection mechanism for otherwise unstable 
systems. 

The role and behavior of resonances during plan- 
etary migration is beyond the sco pe of this cha pter. 
However, the re ader is referred to iPeale (Il976h and 
Malhotral (Il998h for Hamiltonian descriptions of reso- 
nances, which can cleanly treat migration. For a re cent 
applications to exoplanets, see dLee & Pealell2002h for 
how two planets can captu re into a reso nance if their 
migration converges and (iChiandl2003h for how two 
planets can excite each other's eccentricities as they 
pass through a resonance while their orbits diverge. 



Fig. 3. — Orbital element changes after many con- 
junctions. The hypothetical system is the same as 
in Figure |2] except both planets start with eccentric- 
ity, with the inner planet's argument of pericenter 45° 
ahead of the outer planet's. The vertical axis on each 
panel has the same scaling as the corresponding panel 
in Figure 121 emphasizing that the semi-major axes ex- 
perience no net drift, but the eccentricities do. 



3.3. Advanced Interactions 

Having surveyed the basic interactions between two 
planets in we now introduce several more ad- 
vanced topics. 

We previously saw that the eccentricities of plan- 
ets outside mean-motion resonance can change on a 
long timescale. We now extend that concept to sys- 
tems with three or more planets, systems in which the 
orbital elements of each planet vary on many differ- 
ent timescales. Resonances between these timescales 
can excite eccentricities to very high values (see 



Moro-Martm et alJ 120071 for an example in an exo- 



planetary system). For the three-planet system Up- 
silon Andromeda {v And) in the Newtonian approx- 
imation, the innermost planet's eccentricity period- 
ically reaches ^ 0.25, compared to ^ 0.06 in the 
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Fig. 4. — Orbital element changes induced by a mean 
motion resonance. The system is the same as in Fig- 
ure [3] except the outer planet starts at a different pe- 
riod, with a ratio of osculating periods of 2.01 : 1. 
Note that successive conjunctions (where the lines in- 
tersect in the top panel) occur at nearly the same longi- 
tude, which causes the period change to grow. The bar 
shows the vertical scale of each corresponding panel in 
Figures |2]and[3] emphasizing that the semi-major axes 
and eccentricities are experiencing a large oscillation. 



absence of the outer p lanet or ^ 0.0 25 in the absence 
of the middle planet (^ Barnes 2008h. Also important 
in determining the qualitative behavior of the secular 
dynamics is extra precession, e.g., that supphed by 
relativ i stic or tidal effects (^12. lb : see lWu & Goldreich 



(l20CI2l) : lMigaszewski & Goldziewskil(l2009bl) . In Fig^ 



ure|5]we plot the long-term behavior of the eccentric- 
ities in the v And system considering (1) only New- 
tonian point masses, and (2) an extra force modeling 
the tidal bulge raised on the inner planet. The behav- 
ior of the outer planets is not much different, but the 
effect of the extra precession on the inner planet is 
to detune its pericenter precession rate away from an 
eccentricity-exciting secular resonance. Thus, the low 
current value of Cb argues an add itional precession is 
active (lAdams & Laughlin 2006 ): both relativity and 
tides probably contribute with roughly equal preces- 
sion rates of w ~ lO^^^s^^ apiece. 

An extension of secular evolution theory can be 
made into the regime of high inclination and eccentric- 
ity, in which they are strongly coupled (Kozai 196^ 
Lidovl[l962 ). In a system with a planet on an initially 
circular orbit, a third body on a distant exterior orbit 
will periodically pump the planetary eccentricity to a 



maximum of 



v/l - (5/3)cos2i, 



(19) 



where i is the mutual inclination which must initially 
be in the range 39.2° - 140.8°. (For a detailed de- 
scription and a de rivation of this behavior, calle d Kozai 
oscillations, see iFabrvckv & Tremaind 120071 ) Note 



that an initially perpendicular orbit (i = 90°) leads 
to an eccentricity of unity. In some systems, rela- 
tivistic precession would suppress this behavior, but in 
others, tidal dissipation would take hold at the eccen- 
tricity maximum and circularize the planet a t a period 
of a few d ays (Fabrycky & Tremaine 200^ IWu et al. ' 



200?!: iNagasawa et al.ll2008l) . An important indicator 
of whether close-in planets have undergone such an 
event is their orbital orientation relative to the stellar 
spin. With transiting planets, this can be measured 
via the spectroscopic Rossiter-McLaughlin effect (see 
^441 and Chapter 4). 



Ano ther secular effect was found by lMigaszewski & Gozdziewski 
(12009a'). who used a numerical Hamiltonian approach 
to classify stationary solutions for two-planet, non- 
coplanar systems. They found that interlocking plan- 
etary orbits could stably exist even if no mean-motion 
resonance protection is active. Instead, each orbit 
nodally precesses because of the other orbit, and 
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Fig. 5. — Eccentricities as a functio n of time for 
the And system. Orbital elements of I Wright et al. 
(l2009h . and edge-on, coplanar orbits were as- 
sumed. We use nii, — 1.27m0 and (mb, rric, 
md)=(0.672,1.92,4.13) mjup. The outer planets, c 
(Pc = 241.33 d) and d (Pd = 1278.1 d), were 
initialized using Jacobian coordinates, and the plot- 
ted eccentricities are also Jacobian. Two integrations 
were performed, one in which the planets were point 
masses and followed the Newtonian equations of mo- 
tion (eq. 151), and one in which f = fr of equations ^ 
and (O to model the tidal bulge of planet b (Pb = 
4.617136 d), with Rp = 1.3i?jup, A:l = 0.34 as- 
sumed. In the point-mass model, eb quasi-periodically 
reaches 0.25 (top curve of top panel). Including 
the tidal model significantly suppresses the induced eb 
(bottom curve of top panel), and including relativistic 
precession as well (not shown) suppresses it still fur- 
ther The evolution of Cc and are almost identical in 
the two cases, so they are plotted only once. 



they do so at the same rate, so the two orbits re- 
main in the same configuration with respect to each 
other as the whole system precesses. However, in 
numerical tests of the full equations of motion, we 
find that a specific example of this class (fig. 17 of 
Migaszewski & Gozdziewskl2009al) becomes chaotic 



within several Myr unless the mass ratios are smaller 
than~ 10^*. A system of Neptune-mass planets could 
thus exist around a solar-like star in this interlocked 
configuration. 

The high eccentricities of exoplanets call for a mix- 
ture of numerical and analytical work to understand 
the nonlinear dynamics. The typical method of ex- 
panding in small eccentricities and inclinations has 
been fruitfully supplemented with nu merical methods 
(e.g., iMichtchenko & Malhotrall2004l) to average over 
short- timesc ale effects like those of Figure [3] Reso- 
nant orbits_ayii£heccenQn£i^^ 
ties (iLee & Pealel2002l;lBeauge et alj2003l;lLed2004 . 
One way exoplanets can be dynamically simpler than 
Solar System planets is that they are often hierarchi- 
cal, i.e., the period ratios can be large. When this is 
the case, an expansion in the semi-major axis ratio be- 
comes a useful analytic method dFord et al.. .2000, and 

a- 

Mean motion resonances may be shared among 
more than two planets with critical angles which are 
extensions of equation ( fTsT i: these are called three- 
body mean motion reso nances. They h ave been shown 
to destabilize asteroids ( Nesvornv & Morbidellill998t). 
the ou ter planets of the Solar System (IMurrav & Holman 



1999|X and potentially exoplanetary systems (IGozdziewski et al 
.2008ah . 

3.4. Chaos 

Nothing ensures that planets emerging from a pro- 
toplanetary disk — in which collisions were frequent — 
will be on nearly-Keplerian, stable orbits. On the con- 
trary, models of planet formation suggest that the or- 
bits of planets should be closely packed together, with 
the timescale of collisions or ejecti ons being com para- 
ble to a system's current age (e.g., 'L askai boOOV In- 
deed, obse rved multiplanet sy stems are often close to 
instability (IBarnes et al ]l2008h . A dynamical system is 
said to be chaotic if trajectories that are initially sepa- 
rated by an infinitesimal amount diverge exponentially. 
This "microscopic" definition is often used as a com- 
putational tool to detect whether and when "macro- 
scopic" events (i.e., ejections and collisions) are hkely 
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to oc cur dLecar et alj 1 1992c iMorbidelli & Froeschle 
19951) . which usually occur on a timescale orders of 
magnitude longer. 

Res onance overlap is a general condition for strong 
chaos ( IChirikovlll979h . Qualitatively, each resonance 
allows the system to explore a finite zone of semi- 
major axes, called the width of the resonance. If two 
resonances overlap, the system has access to a wider 
swath of semi-major axes. This region might ev e n con 



nect to infinity, allowing an ejection. IWisdom (Il980l) 



found that resonance overlap could account for the 
orbits of massless particles becoming chaotic near a 
planet on a circular orbit. For a planet with a semi- 
major axis of 1 and a mass (relative to the star) of /i, 
he found that for semi-major axes within |Aa|/a < 
1.3 /i^/^ of the planet, a particle's orbit is chaotic: its 
fate is either a planetary collision or an ejection. The 
mass scaling of this limit is found by comparing the 
widths of the resonances to their spacing. 

19891) plementary 



Duncan et al. 



understanding of the chaotic zone. Suppose that two 
planets follow Keplerian orbits except at their mutual 
conjunction, where they give each other a kick. The 
kick leads to a change in eccentricity which is second- 
order in the orbital separation, and a change in orbital 
period which is even higher order If the kick in or- 
bital period is strong enough such that the next con- 
junction occurs more than half an orbital period earlier 
or later than it would have without the kick, the re- 
sult of successive kicks is uncorrected and the system 
executes a random walk, eventually leading to orbital 
instability. This argument leads to nearly the same 
scaling law of the chaotic region: |Aa| < 1.24 /i^/^. 
This c oncept can be re a dily ex tended to N-planet sys- 



tems. [Chambers et al.l (Il99d) applied it to terrestrial 
planet formation (rather small fj, values), found that 
there is no critical separation beyond which stability 
is assured, and mapped out numerical timescales to in- 
stability as a function of orbital separation. Zhou et aTl 
(120070 extended these arguments to the mass ratios of 
most known exoplanets (/i ^ lO"'^), finding that kicks 
at conjunction can explain the instability of planetary 
systems, but that empirical corrections are needed for 
quantitative agreement. 

At the relatively larger masses of most known exo- 
planets (fj, ~ 10"'^), both these methods start breaking 
down: (1) their chaotic zones extend to low p among 
the (p + 1) : p resonances, so the peculiarities of in- 
dividual strong resonances need to be taken into ac- 
count in the resonance overlap picture, and (2) the 



physics of interactions between planets is not nearly 
as localized to conjunctions, and this derivation, based 
on the limit /i — > 0, breaks down. Interestingly, the 
method of resonance overlap no longer breaks down 
when the concept is applied to even higher mass ratios: 
Mardling (2008) has shown that for comparable-mass 
triple stars, n : 1 resonances are the only ones rele- 
vant for stability. This method has been used to find 
the instability bounda ry for planets orbitin g one star of 
a binary star system dMudrvk & Wull2006l) . 

A more general framework for the stability of three- 
body systems is called Hill stability, in which the or- 
bits of two planets can never cross if a particular in 



equali ty is satisfied (Marchal & Bozislll982tlGladman 
.19931: IVeras & Armi tage 2004). There is no known 
sufficient condition for inst ability, unfortunately, but 
iBarnes & Greenbergl20 07b\ave shown that Hill's sta- 
bility criterion and practical stabihty are rather close, 
numerically. 

Chaotic trajectories are particularly important dur- 
ing planet formation, in which a given surface density 
of material is converted into fewer numbers of larger 
bodies. Chaotic zones around a given body scale as 
only a shallow function of its mass, so by putting the 
same surface density into more widely spaced bodies, 
the system becomes more stable. This concept under- 
lies the derivation of the isolation mass — the mass at 
which a protoplanet has cleared its feeding zone and 
no longer accretes — laid out in §2.4 of chapter 13 of 
this volume. 

One popular application of chaotic trajectories 
is the origin of eccentricities. Exoplanets tend to 
have larger eccentricities than the Solar System plan- 
ets. Perhaps they typically start in systems of sev- 
eral planets of comparable mass, which perturb each 
other into crossing orbi ts, resulting in ejection or ac- 
cretion (R asio & Ford 1 19961: IChatteriee eTall 12008 



Juric & Tremainell2008h . If indeed this mechanism ex- 



plains the observed distribution, as discussed in §3.4 of 
chapter 9, it implies that most systems of giant planets 
will self-destruct and ma ny free-floating planets exist. 
See Ford & Rasiol ( 20081 ) for a thorough review of the 
work on this hypothesis and references to competing 
hypotheses for the origin of eccentricities. 

3.5. Stable orbits near planets and in habitable 

zones 



We have seen in that close to a planet, orbits 
become unstable. However, there are particular con- 
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figurations that allow for long-term stability. Although 
none of the following configurations have been ob- 
served yet in exoplanetary systems, it is worthwhile to 
discuss here thek existence (their potential detectabil- 
ity is discussed in [jSTTJ- For each type of configura- 
tion, we consider stable orbits of test particles in orbits 
similar to a planet, then extend the notion to a pair of 
planets. 

Particles could orbit the planet on the orbit of a 
satellite, within its Hill sphere. Consider a parti- 
cle between the star and a planet on a circular orbit, 
in a frame that corotates with the planet. The Hill 
sphere is the region within which the gravity from 
the planet {Gmpr^^, where Vp is the distance to the 
planet) is comparable to the tidal gravity from the star 
(2Gm*rpa~^, in the radial direction). Accounting for 
a differential centrifugal force {Grai,rpa~^) yields the 
customary definition for the Hill radius: 



or 60° behind the planet in its orbit; the se are the sta 



1/3 



(20) 



Bodies orbiting within so me fraction of the Hill radius 
orbit stably as satellites (IDomingos et al.l 120061) . The 
same types of orbits exist even if the satellite is mas- 
sive, even up to the mass equal to that of the planet. In 
this case, the planetary mass nip used to define the Hill 
radius ( eg. EOl) would be th e sum of the two bodies' 
masses (iHenon & Petidll986l) . 

Orbits can also stay near the planet if they are in a 
1:1 resonance, of which there are three kinds. 

The first is an extension of satellite orbits to out- 
side the Hill sphere, cal led quasi-satellite orbits (e.g., 
Shen & Tremainell2008h . From the planet's perspec- 
tive, the particle would look like a very distant satel- 
lite, orbiting the planet in a direction retrograde to their 
common motion around the star From an astrocentric 
perspective, the particle's orbital phase is similar to the 
planet's orbit, but its eccentricity is different, which 
carries it between the planet and the star at periastron 
and to the far side of the planet from the star at apas- 
tron. If that body is endowed with mass as a second 
planet, then both orbits will evolve due to their mu- 
tual perturbation; the eccentric ity can be passed back 
and f ort h between the planets (ILaughlin & ChambersI 



2002h . iHadjidemetriou et al. (l2009l) have explored 



the relationship between this type of 1 : 1 resonance to 
satellite- type mutual orbits. 

The second type of 1:1 resonance is a Trojan orbit, 
named after the asteroids that inhabit this resonance 
with Jupiter Their average position is 60° ahead of 



ble Lagrange points labeled L4 and L5 (iDvorak et al. 
2004). However, the orbits may also wander stably 
around those points, tracing out a shape of a tadpole, so 
they are sometimes called tadpole orbits. Second plan- 
ets of any mass relative to the more massive planet can 
exist in these poin ts, but n = (mi +m-2]lmi, m ust not 
exceed ^ 0.038 (Laug hlin & Chambersll2002l) . or the 



system will be unstable. In this case, both planets will 
trace out tadpole shapes in the frame rotating with the 
long-term mean angular velocity, with the size of the 
tadpole inversely proportional to each planet's mass. 

The width of the stable tadpole region scales as 
/i^/'^, which is steeper than the Hill sphere's scaling 
^^/^ (eq. II20I ). so for low planetary masses there is a 
region between them. In this region lies the third type 
of 1:1 resonance: the horseshoe orbits. Such orbits 
trace out a horseshoe shape in the frame rotating with 
the orbit of the massive planet, which encompass both 
L4 and L5. As with the other resonances, any relative 
mass of two planets may be in this resonance. In fact, 
the Solar System furnishes an example of a 4 : 1 
mass ratio in this reso nance: th e Saturnian satellites 
Janus and Epimetheus (iDermott & Murray 1981.) . 

In systems with known giant planets, the orbital sta- 
bility of hypothetical terrestrial planets has been stud- 
ied extensively, particularly in the region that allows 
for habitable climates (see chapter 16 regarding hab- 
itable zones). An appropriate approximation is that 
the terrestrial planet is too small to affect the orbits 
of the giants: they are treated as test particles (e.g.. 



Sandor et al.ll2007h . Habitable planets might also re- 



side i n the dynamical nich es described above: satel 



lites (IWilliams et al.lll997h or Trojans (iDvorak et al. 



2004 "). In systems already known to have multiple gi- 



ant planets, the zones of stability can be quite compli- 
cated, and numerical integrations are indispensable. 

Terrestrial planets that avoid ejection may still 
be s ubject to osci llating orbital a nd spin properties 
(see Menou & Tab achnik 20031 and Laskar & Robutell 



19931 respectively), and these oscillations may cause 
climate changes (a_ much-enhanced form of Mi- 



lankovitch cycles; iHavs et ak 1976 ). The planet may 
not need to stay in the habitable zone for all of its orbit, 
depending how well t he atmosphere buffers seas onal 
temperature changes (IWilhams & PoUardl l2002r) . so 
the upper limit on a habitable eccentricity is a function 
of planetary properties. It has been argued that moons 
of planets with eccentricities as high as 0.69 (16 Cyg 
B b) might still be habitable. Therefore, ejection may 
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be the only dynamical effect that will spoil a habitable 
world. 

4. HIGHLIGHTS: DYNAMICS IN NATURE 
4.1. GJ 876 and Mean Motion Resonances 

The only exoplanet system hosted by a main se- 
quence star for which non-Kepleri an motion has been 
conclusively detected is GJ 876 dMarcvet all 1200 ll 



Laughlin & Chambersl2001 : Rivera & Lissauerl2001 ). 
The M-dwarf primary hosts three planets, whose prop- 
erties are listed in Table |2] 

These numbers come from a self-consistent Newto- 
nian fit of the radial velocity data (Fig. |6]l, in which 
a coplanar configuration is assumed and the best- 
fitting common sky -inclination is found to be i = 50° 
(IRivera et al.ll2005h . For typical planets discovered by 
radial velocity, which are not found to transit, the data 
are only sensitive to the gravitational influence of the 
planet in the line-of-sight direction, which yields the 
quantity TOpSini, not rrip or i independently. In this 
system, the periastron of each planet rotates at a rate 
proportional to the mass of the other planet, and that 
rate is well determined by the data. Therefore nip is 
independently measured, breaking the irip sin i degen- 
eracy to yield i. 

However, in this system the dynamically-derived 
inclination is in tension with the astrometric orbit 
of planet b dBenedict et aL 2002h . Nevertheless, 



Bean & Seifahrt (l2009h have fit both datasets simulta- 
neously, and they even determined a mutual inclination 
of 4>f,c — 5.0° ^2 3°' close but marginally inconsistent 
with coplanar 

Almost no dynamical constraint can currently be 
given on the orbital orientation of planet d, from ei- 
ther stability considerations or fits to radial velocity 
or astrometry data. Its short period means it is only 
weakly coupled to the outer planets on the timescale 
of the data. 

The dynamically interesting aspect of the system 
is that the outer two planets are deeply engaged in a 
2:1 resonance, with small-amplitude libration of both 
critical arguments 0i = Ac — 2Ah — Wc and 02 = 
Ac — 2A{, — Wb about 0° (compai-e eq. ifTsl ). As a 
consequence, Azu = zuc — Wb = 62 — Oi, also li- 
brates around 0°. The 3-planet, Newtonian, copla- 
nar model with i ~ 50°, for which system parame- 
ters were quoted above, yields |0i|max = 5.4° ± 0.9°, 



jRivera et al.ll2005h . 

Apart from GJ 876, about 7 other planetary sys- 
tems have been shown to be in resonance (Table 143] ). 
This is usually accomplished, for planets of lower 
quality data or fewer dynamical times, by notic- 
ing the sy stem would be unstable if not for the 
resonance ( Vogt et al. 2005a : Correia et alJ 2005at 



Lee et all 12006': 'Fabrvckv & Murrav-Clav"2010'; for 
an overview, see Barnes & Greenberg 2007b and Ta- 



\02 



19.5°±3.8°,and|ATi7|. 



19.4° ±4.3° 



ble l4.1l ). This logic is needed because (1) not enough 
orbital periods have been observed for them to be mea- 
sured to a precision such that non-resonant systems are 
inconsistent with the radial velocity data, and (2) usu- 
ally only a small fraction of the resonance libration 
cycle, or precession cycle, is completed during the 
observation span, so the orbits appear Keplerian. 

4.2. Planets around Pulsar 1257H-12 

The discovery of a planetary system around a 
pulsar (see chapter 8) preceded the first solid ra- 
dial velocity detections around main sequence stars 
dWolszczan & F rail 1992). The discovery technique 
was to infer the gravitational influence of the planets 
through tracking the line-of-sight motion of the host 
star, as in the radial velocity technique. In this case, 
however, the light time effect delays pulse profiles. 
Nevertheless, with only Keplerian motion detected, 
the typical nip sin i degeneracy held sway. 

The theorists set to work quickly, showing that (a) 
the perturbations could be detected with more data, (b) 
they would confirm the planetary nature of the timing 
residuals, and (c) the amplitude and the character of 
the detected perturbations would determine the masses 
of the objects (Rasio et al. 1992: Malhotra et al. l992i) . 
The final point is analogous to how the nip sin i degen- 
eracy is broken for the GJ 876 system (i )4.1l l. However, 
here the main o bservable effect was shown to be pe- 
riod variations (lMalhotralll993bt |Pealelll993b . which 
build up to a sinusoidally varying phase shift with a 
period of 2n/[nB — (3/2)?ic] ~ 5.5yr (see Fig. [TJ. 
The same type of near-resonant behavior is famously 
active in the Solar System between Jupiter and Saturn, 
which has a large effect on their orbital phases known 
classically as the Great Inequality. 

The observers answered the challenge (see chapter 
8), not only detecting pert urbations to c onstrain the 
masses of the planets (Wo lszczanlll994l) . but show- 
ing their orbits are consistent with c oplanar to within 
~ 13° (.Konacki & Wolszcza"nll2003b . The data and the 
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Table 2 
The GJ876 System 



planet 




P[days] 


e 


b 


2.530 ± O.OOSmjup 


60.83(2) 


0.0338 ±0.0025 


c 


0.790 ± 0.006mjup 


30.46(2) 


0.2632 ±0.0013 


d 


7.53 ± 0.70TOEarth 


1.93774(6) 


O.O(assumed) 



No te. — Three-planet Newtonian solution from dRivera et al. 



20051), which optimally fits the radial velocities (fig. |6ll. Uses 



Gm^ = 0.32Gtoq. 



GJ876 ^-planet i = 50° Fit 
Keck 



n I I I I I I 

RMS =, 4.23 m s * 
Vy2 = 1.061 




n 1 1 r 



T 1 1 r 



-J I I L_ 



J I L 



1000 2000 

t (JD - 2450000) 



3000 



Fig. 6 . — Radial velocity data and best three-planet Newtonian fit to the planetary system GJ876. From lRivera et al, 
(l2005h . 
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Fig. 7. — Evolution of the orbital elements of planets B and C around Pulsar 1257+12. The latest orbital elements 
jKonacki & W ol szczan 200 3). taking the planets to be coplanar, are used, and planet A is not included. Following 
Rasio et al.l ( 1992h . the elements displayed are argument of pericenter (cj), change in eccentricity (Ae), and change 
in period (AP) of the Jacobian coordinates of each planet. Error bars of each planet's P from the discovery paper 
(fWolszczan & Frail, 1992) are indicated in the bottom boxes, showing period changes were nearly detectable (but the 
changes in w a nd e were far from detectable). With equally precise data spread over three years, the period variations 
were detected ( Wolszczarilll994 ). 
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PSR B 1257+ 12, Arecibo, 430 MHz 
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Fig. 8. — Time-of-arrival (TOA) of pulses for the pulsar that hosts a three planet system, (a) The daily-averaged 
data with the best triple-Keplerian model as a solid line; (b) the residuals of the Keplerian fit and the difference 
between the best Newtonian and Keplerian models as a solid line, (c) the residuals from the Newtonian model. From 
Konacki & WolszczanI (l2003h . 
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Table 3 
Resonant Systems 



system 


pla^nets 


resonance 




N inner orbits observed 


7Tt^ ^v, 777,^ iTJI T.,^ / sin ?1 
"^p,in; "^p,out L"^Jup/ ^^^^ f'\ 


Ref. 


HD 45364 


b-c 


3 


2 


227 


7 


0.19, 0.66 


A 


GJ876 


c-d 


2 


1 


30 


- 80 


0.56, 1.94 


B 


HD 82943 


b-c 


2 


1 


441.2 


8 


2.01, 1.75 


C 


HD 128311 


b-c 


2 


1 


448.6 


6 


2.18,3.21 


D 


HD 73526 


b-c 


2 


1 


188.3 


13 


2.9, 2.5 


E 


HD 160691 = /i Arae 


d-b 


2 


1 


310.5 


9 


0.52, 1.68 


F 


HD 60532 


b-c 


3 


1 


201 


4.5 


3.15,7.46 


G 


HD 202206 


b-c 


5 


1 


255.87 


8 


17.4, 2.4 


H 


HD 108874 


b-c 


4 


1 


395 


6 


1.36, 1.02 


I 


55 Cnc 


b-c 


near 3 : 1 


14.65 


- 300 


0.82, 0.17 


J 



Note. — References: A: Surrou n ded by cha os /Correia et al."2009t^ B: Ki ng of the r esonant planet s (i^4.11l. C: iMavor et al 



(12004; iFerr az-Mello et al.l (l2005l ) : iLee et alj 12006): Gozdziewski & Konacki (.20061) : iBeauge etal] (l2008bl). D : IVogtetal 



(2005b); Goz dziewski & Konacki 



(l2006h 



E: 



Tinnev et al. (l2006l). F: also known as H D 160691: lPepe et all (l2007h : IShort et al 



G: iDesort et alJ ("2008*): Laskai- & Correial (I2009I) . H: ICorreia et alj (l2005bl) . I: not necessarily resonant dVogt et al 



\). J: according to lFischer et al.. (,2008<) . it is actually just outside this resonance. 



timing model are shown in Figure |8] 

4.3. Secular Apsidal Alignment 

Eccentricities oscillate and apses precess due to 
their planets' secular interaction. If two planets have 
apses precessing at the same rate, on average, they 
are said to be in apsidal lock, with a critical angle 
Aw = zui — ZU2- This angle can librate around ei- 
ther 0° (apses aligned) or 180° (apses anti-aligned), 
depending on the masses and initial orbital elements, 
with the restoring torque supplied by the secular terms. 
(Mean-motion re sonance terms can a lso result in libra- 
tion of Azu, e.g-. lBeauge et al.ll2003l and ^ KT\ but this 
phenomenon is not our current focus.) If one planet 
periodically reaches e = 0, we follow Ford (2008) in 
calling such a system borderline: for such systems, 
Avj is on the border between librating and circulat- 
ing. We note that there is little dynamical signifi- 
cance to this border, and systems on either side of it 
remain close to each other in phase space. In polar co- 
ordinates, (e cos An7, e sin Aw), the difference is just 
whether the trajectory contains the origin or not. 

Recently, there have been several attempts to use 
the libration amplitude, or the proximity to the border- 



line state, to shed light on earlier epochs in multiplanet 
systems. The idea is that if all the planets of the system 
start out in circular orbits, but one is forced to an ec- 
centric orbit, the system will have a small libration am- 
plitude if this forcing is much slow er than the secular 
timescale ( Chiang & Murrayl2002 ), or will be left near 
the borderline state if this forcing is much faster than 
the secular timescale (Malhotra 2002). The agent im- 
parting the initial eccentricity might be the protoplane- 
tary disk or an addi tional planet, which is subsequently 
ejected. iFord et al. (2005) presented the borderline 



behavior of the more massive planets in v And (see 
Fig. as evidence of the latter Later work showed 
that, while borderline behavior i s surprisingly com- 
mon a mong multiplanet systems (iBarnes & GreenbergI 
2006 1). it is perhaps even too common for simple mod- 
els o f scattering among planets to be the explana- 
tion ( Barnes & GreenbergI 2007ah . The final answer 



awaits a rigorous statistical comparison between the 
secul ar behavior resulting fr om scattering simulations 
(e.g., Chatteriee et al. 2008h and the secu lar behavior 
inferred from data (e.g. JVeras et al ]l2009h . 
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4.4. Kozai oscillations 

One of the first very eccentric exoplanets, 16 Cyg 
Bb, caused excitement because its orbit is quite un- 
like the nearly circular orbits of the giant planets 
of the Solar System, and it was unclear how a gi- 
ant planet could form on an eccentric orbit. It also 
suggested that the low eccentricities of the previous 
discoveries were perhaps not primordial, but tidally 
damped. This particular system has a distant compan- 
ion star, whi ch could cause Kozai eccentricity oscil- 
lations (ij2.3. lHolman et al.lll997l:lMazeh et al.lll997l) . 
resulting in the high eccentricity. 



As pointed out by iTakeda et alj ( 12008!) . the four 
planets with the highest eccentricities (e > 0.8), and 
^ 50% of the 18 planets with e > 0.6, have confirmed 
stellar companions; this correlation is statistically sig- 
nificant, considering that onl y 16% of exoplanet 



hosts have stellar companions dMugrauer & Neuhauser 
12009 ) due to survey biases. These numbers may be 
interpreted as a statistical detection of Kozai oscilla- 
tions, which have much too long a timescale to be 
directly observed in the current data of any particu- 
lar system. However, there is additional evidence that 
one of these four systems, HD 80606b with e=0.93, 
has indeed experienced Kozai cycles. It was shown 
that a natura l prediction of the K ozai scenario (as 
described by Wu & Murray 2003 ^ is that the stellar 
spin is currently misaligned from the planetary orbit 
( Fabrvckv & Tremainell2007 ). Th is prediction of mis- 
alignment was recently verified dMoutou et al. 1 120091: 
Iwinn et al. 2009) . 

5. FUTURE PROSPECTS 
5.1. Searching for small planets 

One of the applications of dynamical calculations, 
which has yet to be realized, is the detection of pre- 
viously unknown small planets via their dynamical ef- 
fect on known planets. This potential is particularly 
ripe for the transit timing method. 

Transiting planets offer a wealth of information 
that is inaccessible for a usual radial-velocity-detected 
planet (see chapter 4). An exciting opportunity is af- 
forded by the extreme phase sensitivity of a transit 
hght curve. With high-quahty data, a timing precision 
of '-^ 10s is achievable, which translates to a phase 
measurement of ~ 0.01° along a 4-day orbit. Only 
very slight perturbations lead to along-track variations 
of that order. In comparison, for a Jupiter-mass planet 



on a 4-day orbit around a solar-mass star, a radial ve- 
locity datum with precision of 3 m s~^ translates to a 
phase measurement of ^1°. 

Though it has potential, no compelling examples of 
trans it time variatio ns have yet been published. How- 
ever, ISteffen & Agol (2005 ) — and many authors fol- 
lowing their example — have published stringent up- 
per limits on hypothetical second pla nets, which are 
requir ed by constant transit times 



Agol & SteffenI 



(120070 have analyzed transits of the first-discovered 
transiting planet HD 209458b using the Hubble Space 
Telescope, searching for other companions. They were 
able to put stringent upper limits on the mass of a hy- 
pothetical second planet, as a function of its period: 
see Figure |9] The sensitivity of transit time mea- 
surements is orders of magnitude b etter within reso- 
nances th an midway between them ( Asol et allEoOSl 



Holman & Murray 2005). This is simply a restatement 
of our earlier identification of resonances as effective 
locations for a and P variations. Thus they provide 
complementary information to radial velocity obser- 
vations, which are not particularly good at discover- 
ing planets in resonance (especially the 2:1 resonance; 
Anglada-Esc ude et al. 2010). but they have sensitivity 
over a wide range of periods. 

Although mid-transit times are the most sensitive 
characteristic to period changes over timescales short 
compared to the observations, transit duration varia- 
tions (TDV) have been recognized as more sensitive to 
variations of much longer timescale ( Miralda-Escude 



2002l:lHevl & Gladmani2007l:IPal & Kocsisl2008h . and 
they may eliminate degeneracies inherent to TTV mea- 
surements (lKipping.,2009) . 

Transit timing measurements are a possible way to 
find objects in qualitatively different orbits than those 
available to the radial velocity method. Trojans are 
hard to pull out of the radial velocity data because 
their orbits have the same harmo nics as the main planet 
( Gozdziewski & Konackil 20061) . However, if the two 
planets have low eccentricity and a large libration am- 
plitude, the signal would be a single sinusoid with an 
amphtude that slowly oscillates (every ~ 10 orbital 
period s), markedly different from the signal of a single 
planet (ILaughlin & Chambersll2002l) . Such a libration 
could be easily seen in transit data dpord & HolmanI 
20071) . Trojan planets that make a perfect equilateral 
triangle with the star would not librate at all. How- 
ever, a combination of transit and radial v elocity data 
could still detect it dPord & Gaudill2006h . which has 
been used to place upper limits of varying sensitivity in 
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Fig. 9. — The 3-ct upper limits on the mass of a hypo- 
thetical second planet (with assumed eccentricity 0.02) 
relative to its host star in the HD 209458 system, based 
on the precise and non-variable transit times of planet b 
(solid, gray curve) and that data combined with the ra- 
dial velocity time-series (solid, black curve). A pertur- 
bation theory calculation (thin, dashed curve) matches 
the transit time constraints between resonances, and 
an analytic expression (large dots) matches within res- 
onances. The most sensitive constraints are within res- 
onances, and extend an order of magnitude lower than 
both upper limits from radial velocity (thick dashed 
hne) and the mass of the Earth (thin dashed horizontal 
Une). Near the period ratio of 1, orbits of small bodies 
are withi n the chaotic zon e (j l3.4l ) and are thus unsta- 
ble. From lAeol & Steffenl(l2007i) . 



25 systems jMadhusudhan & Winnll2009h. Moons of 



giant planets could also be searched for ( Simon et al. 



2007 ; Kippind2009 ) by transit timing, but the conceiv- 
able limits are well above the mass of Earth's moon 
or of the moons of giant planets in the solar system. 
The masses of m oons may be hmited b y their forma- 
tion mechanism dCanup & WardI 120061) or by the re- 
quirement that tidal evolution as not destro yed them 
(iBarnes & O'Brienll2002l:[Cassidv et al.1l2009i) . Never- 
theless, many of the giant exoplanets are in the habit- 
able zone of their stars, so the first such planet to transit 
will be carefull y scrut inized for both habitable moons 
(fwilliams et all 1997 ) and habitable trojan compan- 
ions jDvorak et alj|2004 . 



5.2. 



System architectures through detecting per- 
turbations 



In this section we shall examine how more observa- 
tions of non-Keperian motion will contribute to our un- 



derstanding of the architectures of planetary systems, 
some aspects of which we currently know very little. 
We shall see that (a) resonant orbits and (b) mutual 
inclinations of planetary systems will likely enjoy fun- 
damental advances in the coming years via detected 
perturbations. 

First, a historical look at the two well-observed non- 
Keplerian systems shows their enormous contribution 
to our understanding of planets. In 1992, detecting 
perturbations for the planets of PSR 1257H-12 was de- 
cisive in demonstrating the orbital and planetary na- 
ture of th e signal, rather than an unforeseen pu lsar os- 
cillation l^ioietaDliiil l^lszczaSliiii). Thus 
the confirmation of the first exoplanets, and the de- 
tection of the only sub-Earth-mass exoplanet known, 
required the modehng of non-Keplerian motion. In 
2001, detecting perturbations for the planetary sys- 
tem GJ876 was important for demonstrating the true 
masses of planets orbiting main-sequenc e stars; they 



are smaller than masses of brown-dwarfs f Marcv et al. 
2001). Only later did true masses become known for 
many planets thanks to transit measurements. This 
history suggests we can expect perturbations to reveal 
new aspects of planetary systems. 

One aspect of system architectures that is ripe for 
observational input is resonant orbits. The resonance 
of GJ876bc has already been mapped out in detail, 
but orbital changes due to resonance could be de- 
tectab le on decade timescales in other systems as 



well dFerraz-Mello et all 120051: ICorreia et alj 12005a , 



2009h . It might seem that 55 Cnc, with 5 planets, 
two of which are close to resonance, would be a 



good candida te. But neither iMarcy et alJ (120021) nor 



Fischer et alJ (120081) found improved fits when the sys- 
tem was modeled with Newtonian equations ^ rather 
than independent Keplerians. Though perturbations 
have not been detected in most systems, see Table 14.1] 
for a summary of the systems that have been plau- 
sibly claimed to be resonant. Determining whether 
this conclusion is robust, on a case-by-case basis, 
and finding the libration amplitude, has not received 
much observational attention. However, the frequency 
of resonances and the expected libration amplitude 
of the critical angles, or their circulation, has been 
a frequent topic of theoretical research. Currently 
observed resonances have the potential to constrain 
nebular condi tions that lead to t h eir formation and/or 
destruction tLee & Peald 120021; iBeauge et alj 12003 



Adams et al 



20081: 



Lecoanet et al. 



2009h. 



perturba- 



tive events in the system's history (fsandor & Kiev 
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120061; iRein & Papaloizoull2009t iLee et alj|2009l). and 
tidal dissipation for close-in inner planets dNovak et al 



rpla 

|2003; Terquem & Papaloizou 2007|). Therefore, with 
some additional high-quality observational input, we 
stand to learn fundamental information about the ar- 
chitecture of resonances, which in turn informs several 
classes of theories. 

Just outside of resonance, an oscillation in the or- 
bital periods and eccentricities is expected, as shown 
in il4.2l Depending on the proximity from resonance, 
it may be small enough to have eluded detection even 
in well-studied systems, e.g., in 55 Cnc. Conversely, 
since no perturbations are detected there, a constraint 
may be placed on the true masses, as resonance widths 
grow with mass ( Malhotra et a l. 1992). Thus the ab- 
sence of detected non-Keplerian motion can yield im- 
portant constraints. 

A second aspect of system architectures that has re- 
ceived rather little attention is mutual inclination. Al- 
though mutual inclination is recognized as a funda- 
mental quantity for planet formation theory (e.g., it 
was a key inspiration for Laplace's nebular hypothe- 
sis), it cannot be measured by radial velocity if planets 
stay on Keplerian orbits. 

Secular precession due to non-coplanarity has been 



detec ted in triple stars by rad ial velocity (e.g..lJha et al 
2000b and by eclipses (e.g., [Torres & Stefanildl2000l) 



and it has led to constraints on their mutual inclination. 
In a non-coplanar two-planet system, each orbit exerts 
secular torques on the other, and they each precess like 
a top. The precession period in the low-eccentricity, 
nearly coplanar case, with Pout ^ ^in is: 

p ^ ^ 

3 Pin (Pin /Pout) ' 

(21) 

where "in" and "out" refer to the inner and outer plan- 

j II — — I 

ets, respectively jMurrav & D ermott 1999, §7.2). This 

period is also the right order-of-magnitude for moder- 
ate mutual inclinations and eccentricities, and it also 
roughly descri bes the periapse precession of such a 
system jLee & Peale.,2003.) . In systems for which this 
period is small, one might hope to observe precession 
directly and probe mutual inclination. The shortest 
precession period due to secular terms among known 
systems is that of GJ876bc, with Psoc > 100 yr Al- 
though this timescale is much longer than a manage- 
able observing program, even a small part of this cy- 
cle could produce observable effects, because over the 
complete cycle the orbit can entirely reorient. So far 



mutual inclination in planetary systems has only been 
measure d in two systems, based on resonant pertur- 
bation s dKonacki & Wolszczanll2003l; iBean & SeifahrtI 



2009) 



Transit-timing variations on orbital timescales may 
be a much faster route to mutual inclinations. The- 
oretical evidence has been building that the pattern 
made by transit times can reveal mutual inclination 
in a single passage of an external, eccen tric planet 
(iBorkovits et al.1 120031; iBakos et al.ll2009bh or the de- 



tailed si gnal of s hort-term interactions (lAgol et al 
.200 5; Nesvorny & B eauge 2010). Therefore, as soon 
as transit time variations are detected, a constraint can 
be put on mutual inclination. 

The transit-timing method may prove to be crucial 
in interpreting aspects of the HAT-P-13 planetary sys- 
tem, for w hich planet b, the inne r planet, tra nsits. Ac 



cording to iBatygin et al.l (l2009h . following iMardling 



(l2007h . planet b should have damped to a calcula- 
ble, non-zero eccentricity due to forcing by the outer 
planet, which is massive and eccentric. The value of 
this forced eccentricity depends on the planet's tidal 
deformability (the Love number k^) and thus its inte- 
rior structure. Thus for the first time, we could have a 
constraint on the mass distribution interior to an exo- 
planet. However, this chain of logic assumes that the 
two plane ts are coplanar, so i t is crucial to establish 
that fact. Bakos et al. I (l2009bh showed in the discov- 
ery paper that ~ 5 s transit time variations are ex- 
pected if the system is coplanar, and even larger vari- 
ations (of a different shape) are expected if the sys- 
tem is non-coplanar Therefore, a tidal bulge may be 
soon inferred, resting on two aspects of non-Keplerian 
dynamics: (1) orbital timescale perturbations for ob- 
servers to determine the mutual inclination, and (2) 
secular timescale and tidal evolution, which causes the 
interior structure of the planet to feed back on its ob- 
servable eccentricity. 

For further reading, the standard reference textbook 
for planetary dynamics is Solar System Dynamics by 



Murray & Dermott (Il999l) . It contains much of the 



material here, with all the required mathematical de- 
tail, but it was written too early to include many re- 
sults on exoplanets. Two recent reviews on the dy- 
namics of exopl anet systems fill in many of the details 
of this chapter: iMichtchenko et al.l ( 120071) treat sec- 
ular and resonant ef fects in two-planet systems, and 



Beauge et alj (l2008ah discuss resonant dynamics (par- 



ticularly of the 2:1 resonance) and orbital fits to such 
systems. 
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The study of non-Keplerian motion in exoplanets 
has played a critical role in understanding their nature 
and histories, from the first known exoplanets (around 
a pulsar, of all places) to systems of planets in dynam- 
ical configurations which shed light on their early his- 
tory. The relatively simple equations hold a wealth of 
complexity, the extent of which is still being mapped 
in parallel with the continued discovery of planets un- 
Uke those of our Solar System. 



Note added in proof: The first haul of candi- 
date multiple-planet systems discovered by the transit 
method has now bee n brought in by the Kepler mission 
JSteffen et al.ll2010l) . which suggests timing measure- 
ments (§§2.5, 5) will be revolutionized, for several 
reasons. First, Kepler is efficient at detecting much 
longer-period transiting planets than has been possible 
by ground-based surveys, and the timing signals scale 



with t he orbital periods of the planets (e.g.. lAgol et al. 



20051). Besides that simple scaling, systems with pe- 
riods > 10 days are little affected by tidal dissipation, 
so pairs of planets can stay in resonance, resulting in 
very large timing signals their whole lifetime. Second, 
in the single-transiting case, various system architec- 
tures can yield nearly indistinguishable timing si gnals 



jFord & Holmanl2007l:lNesvornv & Morbidellil2008i) : 
but in systems of two planets that both transit and show 
timing variations, the system pa rameters (including 



masses; 



Holman & Murravll2005h can be much more 



easily inferred from the data. Finally, ISteffen et al. 



(12010 ") showed that their 5 systems with multiple tran- 
sits will have timing variations readily apparent in the 
full Kepler dataset. They also showed that for every 
system with multiple planets that transit, due to view- 
ing geometry there must be 2-20 or more systems like 
it, each with only one planet displaying transits. These 
two facts together imply that Kepler will discover tens 
if not hundreds of planets on detectably non-Keplerian 
orbits ! 
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